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ABSTRACT 



CZ3 



Context. Weak flares are supposed to be an important heating agent of the outer layers of stellar atmospheres. However, due to instrumental 
limitations, only large X-ray flares have been studied in detail until now. 

Aims. We used an XMM-Newton observation of the very active BY-Dra type binary star CC Eri in order to investigate the properties of two 
^ ' flares that are weaker than those typically studied in the literature. 

Methods. We performed time-resolved spectroscopy of the data taken with the EPIC-PN CCD camera. A multi-temperature model was used 
■ to fit the spectra. We inferred the size of the flaring loops using the density-temperature diagram. The loop scaling laws were applied for 
[ deriving physical parameters of the flaring plasma. We also estimated the number of loops involved in the observed flares. 
' i . Results. A large X-ray variability was found. Spectral analysis showed that all the regions in the light curve, including the flare segments, 
' are well-described by a 3-T model with variable emission measures but, surprisingly, with constant temperatures (values of 3, 10 and 22 
CZ^ MK). The analysed flares lasted ~ 3.4 and 7.1 ks, with flux increases of factors 1.5 - 1.9. They occurred in arcades made of a few tens 
of similar coronal loops. The size of the flaring loops is much smaller than the distance between the stellar surfaces in the binary system, 
and even smaller than the radius of each of the stars. The obtained results are consistent with the following ideas: (i) the whole X-ray light 



o 



^» , curve of CC Eri could be the result of a superposition of multiple low-energy flares, and (ii) stellar flares can be scaled-up versions of solar flares. 



X 
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1. Introduction 

Solar-like stars (main sequence stars with spectral types from 
F to early M) have radiative cores and convective outer en- 
velopes. Convection, together with differentia l rotation, gen- 
erates a magnetic dynamo (see IParkei , 1975 , and references 
ther ein) that is responsible of the formation of the corona 
(see Favata & Micelal l2003t land Gtidell 12004 for two exten- 
sive reviews about stellar coronal astronomy). Magnetic ac- 
tivity similar to that observed on the Sun is typically detected 
on these stars, showing variability through all the electromag- 
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netic spectrum. The activity level of a star is frequently mea- 
sured in terms of its coronal X-ray luminosity L x . Using data 
collected by t he ROSAT satellite, ISchmitt et all d 1995b and 



Schmitt (1997) found that solar-like stars present X-ray lumi- 



nosities (integrated over the 0. 1-2.4 keV energy band) in the 
range 25.5 < log Lx(erg s -1 ) < 29.5. This luminosity is cor - 



related with the stellar rotation rate (Pall avicini et all 1 1 98 1 ). 



IWalterl |l982) noticed that a single power law dependence be- 
tween Lx/Lbd and the angular velocity was unable to repro- 
duce all the observed data, and proposed to replace it by ei- 
ther a broken power law or an exponential relationship. This 
result later led to the c oncept of saturation of stellar activ - 
ity at high rotation rates dVilhul Il984t IVilhu & Walteri Il987l) . 
Rotational velocity, and therefore activity, decreases with age 
because of the an gular momentum losses due to the magne - 
tised stellar wind (IRengaraianl Il984t IPace & Pasquinii 120041) . 
However, short-period binaries can maintain high rotation rates 
since they are synchronized by tidal coupling. 
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Flares are the most extreme evidence of magnetic activ- 
ity in stellar atmospheres. Frequent flaring is found on late 
K and M dwarfs in the solar neighbourhood (the so-called 
UV Ceti-type stars - see their general properties in Pettersen, 
1991). Similarities between solar flares and those observed on 
UV Ceti-type stars suggest that they are produced by the same 
basic physical mechanisms. Flares are supposed to be the result 
of the energy relea se from magnetic field r econnection in the 
lower corona (e.g. iKopp & Pneuman , Il976l) . In consequence, 
electrons and ions are accelerated and gyrate downward along 
the magnetic field lines, producing synchrotron radio emission. 
Bremsstrahlung radiation is emitted in hard X-rays (> 20 keV) 
when these ionized beams collide with the denser material of 
the chromosphere. At the same time, the gas in the affected 
chromospheric region is heated (optical and UV radiation is 
then emitted) and evaporated. Thus, the density and tempera- 
ture of the newly formed coronal loops increase, emitting in 
soft X-ray and extreme UV wavele ngths. Recent evidence in 
favo ur of this scenario was given by lMitra-Kraev et alJ ([2005) 
and lSmith etall d2005b . 

A recurrent idea in the literature is that flares are the 
main heating agen t of the outer stellar atmospheres (e.g. 
Audard et al. , 2000l) . so that the observed "quiescent" emission 
would be the result of a superposition of multiple small flares 
(named "nano-flares"). The distribution law that gives the num- 
ber of flares dN within an energy-interval [E, E + dE] allows us 
to estimate the energy budget of the corona. I n particular, so 



lar flares are distributed following a power law dDatlowe et al 



1974: ILin et all 1 1984: Dennis, 1985), that is 



dE 



(1) 



where k\ is a constant and a is the power-law index. For a > 2, 
an extrapolation of Eq. (Q} to flare energies below the detection 
threshold would be sufficient to account for the luminosity of 
the quiescent corona. For this reason , it is crucial to investigate 
the validity of such an extrapolation ([Hudson!, 1 199 II) . as well as 
the value of a on magnetically active stars. For the Sun, values 



of a between 1.5 a nd 2.6 were reported (jCrosbv et al. , 1993 



Porter et al., 1995; Krucker & Benzl 



1998T 



However, a later 



study carried out by Aschwanden et alJ ( 2000) suggests the in 



sufficiency of nano-fl ares to heat the solar cor ona. On the other 
hand, for active stars. ICollura et alJd 19881) and lOsten & Brown 
(1999) ? ound a » 1.5 - 1.6, while other authors have re- 
ported a > 2 (I Audard et allll999ll2000l ; lKashvap et all 12002 : 



Giidel et al. , 2003b . Therefore, it is still unclear if Eq. ([]} can be 
extrapolated from large observable flares towards the weakest 
ones in order to support the flare-heating hypothesis. 

While flare stars have been studied in the optical during 
more than half a century, the first sizeable sample of X-ray 
flares was o nly compiled after the launch of the Ein stein 



(European Photon Imaging Cameras) detectors - on-board the 
XMM-Newton satellite - also enable the detection and anal- 
ysis of smaller flares. UV Ceti-type flare stars are specially 
indicated for this purpose because of its proximity. For all 
these reasons, we decided to carry out the present study of the 
X-ray flares detected on the UV Ceti-type star CC Eri using the 
XMM-Newton satellite. 

CC Eri (HD 16157) is a spectroscopic binary star (BY 
Dra-type) located in the immediate solar nei ghbourhood, at a 
distan ce of 11.51 + 0. 1 1 pc (f rom Hipparcos. IPerrvman et all 
1997 ). This is a SB2 system ( Strassmeier et all " 1993b which 
consists of a K7.5 Ve primary and a M3.5Ve seconda 
2000), with mass ratio 



(Amado et al 



Amadoetal 



2 ((Evans, 

20001) . The photometric period - 1 56 day s 



ndary 



suits to be equal to that of the orbital mo tion ( Evans , 19591 : 
iBopp & Evanlll973l:lBopp & FekelLll977b . Thus, the synchro- 
nization due to the tidal lock makes the primary component 
to be one of the fastest rotating late K dwarfs in the solar 
vicinity. Using kinematic al criteria, the age of th e system was 



estima ted t o be 9.16 Gyr dDemircan et al. , 2006). Buskoetal 



d 1977b and lAmadoeTal] (2000) found that the chromospheric 
emission of CC Eri varies in antiphase with its optical contin- 
uum, suggesting the presence of active emission regions asso- 
ciated with starspots. Besides, i ts quiescent radio emission is 



2' 



olarized at the 10 - 20 % level dOsten etall 12002; Sleeetal 



2004), indicating large-scale ordering in the stellar magnetic 
field. CC Eri pr esents a strong flare activity over a wide 



......... . . "tiies di 

19881: iBvrne et all 1 19921 IGiidelL 1 19921: |Pan & Jordanl 1 1995 



range of energies (Busko & Torres, 1976, 1978; Caillault et al 



_ IGiidelL . 

Ama do et all |2000|; bsten et all l2002t Isiee et all l2004h . First 



X-ray detections of CC Eri were done with HEAOl , showing 
log L x(erg s -1 ) « 29.26 in the 2-20keV band dTsikoudil 
119821) . Einstein, which m easured log Lx (erg s ) ~ 29.51 in 
the 0. 15 - 4.5 keV range dCaiHaultL Il982l) . and EXOSAT, that 
obs erved log Lx(erg s~') « 29.62 in the 0.04 - 2 keV inter- 



val dPallavicini et allll988l) . However, since none of these ob- 



servations lasted very long, little information on the temporal 
and spectral variatio n of the source in X-rays was available. 
Pan & Jordanl d 1995b observed and analysed for the first time 
an X-ray flare on CC Eri, using ROSAT observations. They 
measured a quiescent luminosity log Lx(erg s _1 ) « 29.40 in 
the 0.17-2 keV band. The flare had an e-folding rise and de- 
cay times of about 1 h (or less) and 2 h, respectively, and the 
emission was enhanced by a factor greater than 2. A 2-T model 
gave an adequate description to all the ROSAT spectra and 
showed the presence of high-temperature plasma (~ 10 MK) 
even during the time-intervals where no flare activity was de- 
tected. This is consistent with the results obtained for the qui- 
escent emission of M-d warfs using data f r om th e EXOSAT ME 



and E instein IPC (e.g. iPallavicini et al. , . 1990; Schmitt et al 



Observatory (Haisch. ll983l) . Later. lPallavicini et al. d 1990b pre- [l990). All the X-ray luminosities measured for CC Eri are sim- 



sented the results of a comprehensive survey of X-ray observa- 
tions of flare stars carried out with the EXOSAT Observatory. 
However, until recent years, and due to instrumental limita- 
tions, only large X-ray flares could be studied in detail. At 
the present time, the great sensitivity, wide energy range, high 
energy resolution, and continuous time coverage of the EPIC 



ilar and place this binary star among the most active ones. This 
is clearly noticed when comparing the quiescent X-ray lumi- 
nosity obtained for CC Eri from the data analysed in this work 
(Lx, 0.1-2.4 kev ~ 3.7 x 10 29 erg s" 1 ) with the cumulative X-ray 
luminosity distribution functions found bv lSchmitt et all d 1995b 
for low-mass stars in the solar neighbourhood. 
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Fig. 1. Source (upper region) and background (lower region) 
light curves of CC Eri as observed with the EPIC-PN detector 
for the energy band 0.5 - 10.0 keV and a temporal binning of 
100 s. The labels FA, FB and FC refer to the detected flares 
while L0, LI and L2 designate different emission levels where 
no clear flare activity was observed. 

In this paper, we present the study of an XMM-Newton ob- 
servation of CC Eri, which, for the first time, reveals details 
about the behaviour of the parameters that characterized the 
plasma during two flares weaker than those typically analysed 
in other active dM stars. Technical information about the ob- 
servation and details of the data analysis is given in Sect. [2] 
In Sect. [3] we describe the light curve and the different kinds 
of variations observed. The time-resolved study of spectral pa- 
rameters is presented in Sect. |4] where we also estimate the size 
of the flaring loops. Finally, in Sect. [5] we discuss and interpret 
the results in the context of solar and stellar flares. 

2. Observations and data analysis 

The observation of CC Eri analysed in this work was per- 
formed with the XMM-Newton satellite on August, 8, 2003 
(PI: H.Kay, ID: 0148790101) during revolution 0671. The 
XMM-Newton satellite owns the most sensitive soft X-ray de- 
tector system presently available: the EPIC instrument. It con- 
sists of three imaging and non-dispersive C CD-based cameras: 



the twin MPS 1 an d MOS 2, and the PN dTurner et aU 12001 



Strii der et all 1200 II) . The wavelength range measured by the 
EPIC detectors allows one to obtain a reliable determination 
of properties of the hottest plasma components. For our timing 
and spectral analysis, we used only the data from the PN-CCD 
camera, that is more sensitive than the MOS detectors. This 
observation was done in the full frame mode with the thick fil- 
ter. The exposure time of the image taken with EPIC-PN was 
36.703 ks, that is, 0.27 times the orbital period of the binary 
star CC Eri. 

We used the event file in the PPS data products, which was 
produced with the standard XMM-Newton Science Analysis 
System (SAS) software, version 5.4.2. The light curve and the 
spectra were obtained with standard tools of SAS. The PN re- 



Table 1. Time-intervals of the main activity levels detected on 
CC Eri (see Fig. [B- 

Activity level Time-interval (ks after the observation starts) 



L0 


0.0- 3.6 


FA 


3.6- 5.6 


FB 


5.6-12.9 


LI 


12.9-17.7 


L2 


17.7-23.7 


FC 


23.7-36.7 



Table 2. Duration of the rising (t r ) and decay (td) phases of 
the flares observed on CC Eri. 



Flare 


tr (s) 


td (s) 


FA 


400 


580 


FB 


570 


2850 


FC 


1100 


5960 



sponses were generated with the SAS RMFGEN and ARFGEN 
tasks. The spectral analysis was done with the X-ray spectral 
fitting package XSPEC V 1 1 . 3 .2 dArnaudl Il996i \20o4) . 

Stand ard selection cri teria were applied for filtering the 



data (see Ehl e et aUl2004l) . We extracted events in the energy 
band between 0.5 - 10.0 keV that triggered only one or two 
detector pixels at the same time (PATTERN < 4). Data below 
0.5 keV were excluded to avoid residual calibration problems 
in the response matrices at soft energies. Nevertheless, since 
flares mainly affect hotter thermal components, the spectral re- 
gion below 0.5 keV is not crucial for our analysis. In fact, we 
checked that no significant variations are found in the spectral 
parameters (see Sect. 14. Il l when the region 0.3 - 0.5 keV is in- 
cluded for fitting the spectra. The EPATPLOT task was used for 
confirming the existence of pile-up affecting the inner region of 
CC Eri. To lose the minimum number of counts as possible, we 
looked for and ignored the smallest region that allowed avoid- 
ing the pile-up effects during the whole observation. The X-ray 
light curve and spectra were therefore obtained with the events 
taken from an annulus with inner radius 12" and outer radius 
44". On the other hand, background photons were extracted 
from a source-free region - a circle with a radius of 54.4"- 
placed at the same CCD as CC Eri. 

3. The light curve 

Fig. Q] shows the source (background - subtracted) and back- 
ground light curves of CC Eri as observed with the EPIC-PN 
detector for a temporal binning of 100 s and the filtering cri- 
teria given in Sect. [2] Note that the background count-rate was 
scaled to the size of the source region. The dead-time correc- 
tion was applied to both the source and background count-rates. 
A significant X-ray variability was found throughout the obser- 
vation. The lowest activity level (L0 or quiescent state here- 
after) was observed at the beginning. It was followed by a 
small flare (FA) where the stellar flux increased by a factor 
1.3. Another two flares (FB and FC) with flux changes of fac- 
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Fig. 2. Cumulative distribution function of the count-rates observed in the regions LO, LI and L2 (solid black line) compared 
to those simulated by supposing a constant source with count-rate equal to the mean value (/? m [Li]) in the given region (dashed 
area). 




Rote (count/s) Rote (count/s) 



Fig. 3. Cumulative distribution function of the observed 
count-rates obtained for LO, LI and L2 altogether (solid black 
line) compared to those simulated by supposing a constant 
source with count-rate equal to the mean value of these three 
regions (dashed area). 



Fig. 4. Cumulative distribution function of the count-rates de- 
tected in the whole observation of CC Eri (solid black line) 
compared to those simulated by supposing a constant source 
with count-rate equal to the quiescent value (dashed area). 



tors 1.5 — 1.9, depending on the selected reference level, were 
also observed. The maxima of these two flares were separated 
by 5 hours. In addition, the stellar emission between these two 
events was also variable - two different activity levels (LI and 
L2) were identified - and higher than that observed during the 
quiescent state. The start and end times of the main selected 
periods are summarized in Table Q] FA, FB and FC radiated 
in the 0.5 - 10.0 keV band a total energy of 0.08, 0.75 and 
1.5 x 10 33 erg, respectively. 

In terms of relative increasing of flux with respect to the 
quiescent level, the strength of the detected flares is some- 
what smaller than the typic al values analyzed o n acti ve low- 
mass stars. For instance, iRobrade & Schmittl d2005l) found 
flux increases of factors 2 - 3 in flares produced by EQ Peg, 
AT Mic, AD Leo and EV Lac. All of them can be consid- 
ered as moderate flares in view of the fl ux increases ( 10 - 30 
times the quiescent state value) that iFavata et al.l (|2000cJ), 



Katsova et alj 020021) . and lGtidel et a l. (2004) observed respec 
tively on EV Lac, EQ Peg and Prox Cen. The detected flares 



are eve n a little bit weaker than that observed by [Pan & Jo rdan 



d 19951) on the same star, CC Eri, which showed a ratio be- 
tween the peak and minimum fluxes larger than a factor of 2. 
However, in absolute terms, the peak X-ray luminosity of the 
analysed flares is large compared to that of the brightest solar 
flares. 

Table [2] lists the duration of the rising and decay phases of 
the observed flares. Since the behaviour of these two phases can 
be well-described by an exponential law, their duration was re- 
spectively estimated as the 1/e rise time (tr) or 1/e decay time 
(td). The values of tr and td were therefore determined from 
a least-squares fit to the corresponding data by an exponen- 
tial function of the form R = A( ) e ( '~' ,mx) / T + Rb- dS e\ where R is the 
count-rate, Rb Rse is the count-rate in the quiescent state (L0 in 
this case), Aq is the amplitude at the flare maximum, t is the 
time, f max is the time at the flare maximum, and t is Tr or - t d . 



I. Crespo-Chacon et al.: XMM-Newton observation of the flare star CC Eri 



5 




20 
Time (ks) 

Fig. 5. Rebinned light curve of CC Eri for the energy band 
0.5- 10.0 keV. 



Note in Fig.Q]that flare FC could have a larger or shorter decay 
time, depending on the value considered for Rbase (i.e., equal 
to that of L0 or, on the contrary, LI), but being always greater 
than 50 minut es. Using data collect ed by the LE experiment 
on EXOSAT. IPailavicini et all dl990h found two different types 
of flares in a sample of 32 M dwarf stars, i.e., impulsive flares 
and long decay flares, s imilar to the ones observed on the Sun 



(Pall avicini et all 119771) . Those in the first group are reminis- 



cent of solar compact flares, showing rise times of a few min- 
utes and decay times of tens of minutes; and those in the second 
group, with decay times of the order of w 1 hour or longer, are 
reminiscent of solar long-duration 2-ribbon flares. Thus, FA 
can be classified as an impulsive flare and FC as a long de- 
cay one. However, the classification of flare FB is unclear since 
its decay time is in the limit ra nge between the two typ es. As 
discussed bv lPallavicinl dl988l) and lPoletto et alj dl988h . these 
morphological differences may indicate real physical differ- 
ences in the energy release process, as i t also appears to happen 



for solar co mpact and 2-ribbon flares dPallavicini et al. . 1 1 977t 



Priest, 198 1 ). In particular, in compact flares energy is probably 
released only during the impulsive phase, whereas in 2-ribbon 
flares a prolonged energy release is apparently required to ex- 
plain their long decay times. Regarding this point, it is inter- 
esting to note the peaks observed in the light curve during the 
decay of flare FB, which are probably related to different mag- 
netic reconnection processes (superposed flare-type events). 

3.1. Searching for short-term variability 

There are three long intervals in the light curve of CC Eri (L0, 
LI and L2) that have been supposed to have a constant emis- 
sion level. In this se ction we make use of the method given by 
Marin o et al. (2000) to test whether small changes within each 
one of these time regions (see Fig.[TJ can be considered further 
short-term stellar variability or, on the contrary, are compatible 
with statistical fluctuations. 

Assuming that during the intervals L0, LI and L2 the 
star has a constant count-rate equal to the mean count-rate 
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Fig. 6. X-ray spectrum of the quiescent state of CC Eri (L0). 
The best fit (2-T APEC model) is also shown with a solid line. 

characteristic of each one of these regions (/? m [Li]), we com- 
puted the net counts expected (c ex p, bin[Li]) in the used tem- 
poral bin (100 s). For every interval, we generated a set of 
simulated data with the N Poisson distribution centered on 
Cexp, bin [Li], which represents possible outcomes from the ob- 
servations if the source had a constant count-rate equal to the 
mean count-rate of the interval. Each set consisted of 1000 of 
such simulations. We then calculated the cumulative distribu- 
tion function (cdf3) of the count-rates observed within the in- 
terval, as well as of each one of these simulations. The set of 
the CDF's obtained for the simulations allowed us to evaluate 
the spread introduced by statistical fluctuations. 

In Fig. |2] we have plotted the CDF of the intervals L0, 
LI and L2. All of them are compared with their correspond- 
ing simulations for a constant source with count-rate equal to 
R m [Li] (dashed area). In all the three cases the CDF of the ob- 
served count-rates is contained within the space occupied by 
the CDF's of the simulations. Therefore, statistical fluctuations 
can account for the observed spread and possible short-term 
variability on CC Eri cannot be distinguished from noise in 
these regions. However, when the analysis is done for the three 
intervals (L0, LI and L2) altogether, significant differences be- 
tween the observations and simulations are found (see Fig.0. 
In fact, the same happens when the analysis is separately done 
for L0 and LI, or LI and L2. The difference between the emis- 
sion from L0, LI and L2 is therefore due to stellar variability 
and cannot be explained only by statistical fluctuations. 

The same type of study can also be used to quantify the 
variations found in the complete light curve of CC Eri (Fig.|4|i. 
Since in this case we were looking for changes relative to 
the quiescent emission, we chose as reference level the mean 
count-rate measured for the quiescent state, i.e. R m [L0]. If the 
duration of the observation were long enough, the CDF could 
be considered as the fraction of time that the star spends with 
a count-rate greater than, or equal to, a given value. The tail of 
the CDF at rate > 3 count/s is produced by the two strongest 



1 The cumulative distribution function (CDF) represents the prob- 
ability of observing a number of counts - in the chosen temporal 
bin - greater than, or equal to, a given value. 
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detected flares (FB and FC). The rest of the CDF not consistent 
with the simulations accounts for the other kind of variability 
observed at different time scales on the light curve of CC Eri. 

In addition, we analyzed FA to confirm that the emission 
within this time-interval was not constant. We obtained that the 
part of the CDF at the highest count-rates of FA is out of the 
simulations. Therefore, FA is probably a flare-event. However, 
the current instrumentation did not allow us to carry out a 
time-resolved study of the different phases of FA due to the 
low strength and shortness of this flare. 



4. Spectral analysis 

In this section we carry out a time-resolved study of coro- 
nal properties (temperature, emission measure and metal abun- 
dance) through the whole observation of CC Eri. 

With this scope, and taking the variability pattern in the 
light curve into account, the total observing time was split in 
several time-intervals with approximately constant count-rate. 
All the intervals were chosen in order to have enough signal to 
perform a reliable spectral analysis. Most of the spectra have 
at least 2000 counts. However, the spectrum extracted for the 
rising phase of both FB and FC was taken with less photons 
(> 1000 counts) for avoiding to mix this region with the flare 
maximum. The rebinned light curve of CC Eri is shown in 
Fig. [5] We obtained the spectra of all these intervals follow- 
ing the filtering criteria given in Sect. [2] Every spectrum was 
binned to provide at least 8 counts per spectral bin. Bad chan- 
nels were always excluded. 

We modeled the spectra by using the Astr ophysical 
Plasma Emission Code (APEC, ISmith et all l2001al) included 
in the XSPEC software. APEC calculates spectral models for 
hot, optically thin plasmas using atomic data s tored in the 



Astrop hysical Plasma Emission Database (APED. [Smith et al 



2001bh . The APED files contain atomic data such as colli- 



sional and radiative rates, recombination cross sections, di- 
electronic recombination rates, and satellite line wavelengths, 
which constitute the relevant information for calculating both 
the continuum and line emission. Interstellar absorption was 
ta ken into account by using the photoelectric cross sections 
of lMorrison & McCammonl d 19831) . also available in XSPEC. 
Despite the generally satisfactory results of using the x 2 min- 
imization technique for obtaining the best fitted model, it runs 
into problems when the number of events is small. For this rea- 
son, our fit procedure was based on the Cash's C statistic mini- 
mization dNousek & Shue , Il989h . which gives better fits in the 
low-count regime (note also that the C statistic is equivalent to 
X 1 in the limit of large number of counts). Unfortunately, for 
the C statistic there exists no method analogous to that of the 
reduced^ 2 value Ofr ed ) with which we can measure the good- 
ness of the fit. We can determine the best parameters by min- 
imizing the function, but we have no criteria for rejecting the 
model. Therefore, for giving an idea of the goodness of our fits 
we show the reduced x 2 value generated by each model (ob- 
tained through C statistic minimization) and its corresponding 
data set. The errors associated with the fitted parameters were 
calculated for a confidence level of 2.706 cr. 



A 2-T model was needed to describe the shape of the 
coronal spectrum of the quiescent state (L0, see Fig. |6). 
Leaving all the parameters as free variables, and taking the 
same global abundance (scaled o n the solar photospheric val- 
ues o f lAnders & Grevessel fl989l) for all the temperature com- 
ponents, we obtained the following results: the hydrogen col- 
umn density (Nu) of the intervening interstellar medium is 
negligible in the studied wavelength range, as expected due 
to the proximity of CC Eri; the global coronal abundance is 
Z/Z G = 0.33^'°; the temperatures are KTi = 0.301+°°}| keV 
and KT2 = 0.927^^ keV; and the emission measures for 
these thermal components are EM\ = 1.69^2 3* X 10 2 cm" 3 
and EM 2 = 1-50+° ^ x 10 52 cirT 3 , that is, EM 2 /EM\ ~ 0.9. 
These values give the best possible fit to the data using the C 
statistic, corresponding to ax 2 d ~ 1-44. 



4.1. Time-resolved study of spectral parameters 

We carried out the spectral analysis of all the time-intervals 
showed in Fig.|5]to examine the influence of the observed flares 
on the temperatur es and emission measures of the plasma. 
Real e et al. (2001) converted X-ray data of a sample of solar 
flares into the same format and framework as stellar X-ray 
data, in the perspective to use them as templates for inter- 
preting stellar flares. They found that synthesized stellar-like 
spectra of solar flares (previously subtracted by the quiescent 
"background" spectrum) are generally well-fitted with a sin- 
gle thermal component (l-T model) at a temperature close to 
that of the maximum of the EM( T) distribution. However, two 
thermal components are sometimes needed to fit the data dur- 
ing the decay, probably due to a rapid variation of the plasma 
temperature within this phase. In fact, they detected deviations 
from the isothermal description for the flares with no signifi- 
cant sustained heating. Consequently, our first approach con- 
sisted on using a 3-T model to fit each observed spectrum, fix- 
ing all the parameters of the first two components at their quies- 
cent values, whereas the temperature and emission measure of 
the third component (which describes the flaring plasma) were 
free to vary. Generally, the statistics did not allow to constrain 
the abundance of the additional APEC component. Therefore, 
it was assumed to be equal to the global abundance obtained 
for the quiescent state in order to avoid unphysical solutions. 
Results showed that this model failed to reproduce the observed 
data under the considered assumptions. A systematic excess at 
high energies appeared in the residuals, and the temperature of 
the hottest component resulted to be very similar to that of the 
second one. This suggested us that the emission measure of the 
second component may be changing (note that we fixed it at its 
quiescent value), and therefore the third component was trying 
to compensate its enhancement providing a temperature lower 
than that required to fit the high-energy region of the spectra. 
For these reasons, we decided to perform a self-consistent anal- 
ysis by fitting the data using an "iterative" procedure, which is 
described below. 

Each spectrum was fitted with a 3-T model, taking the 
same global abundance for all the thermal components and fix- 
ing it at its quiescent value. The first iteration was done leaving 
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Fig. 7. Time-evolution of the coronal temperatures obtained with the 3-T model. In the KT$ panel, a dashed line marks the 
time-segments for which this temperature could not be constrained (see text). Note that the time-resolution is the same as that 
shown in Fig. [5] 
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Fig. 8. Coronal temperatures (obtained with the 3-T model) vs. count-rate. In the KT3 panel, a dashed line marks the 
time-segments for which this temperature could not be constrained (see text). Note that no significant changes in the temper- 
ature values are observed even at flare maxima (points with the highest count-rates). 
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Fig. 9. Distribution of the three temperatures that characterize the corona of CC Eri during the XMM-Newton observation. The 
corresponding Gaussian fits are also plotted. 



all the temperatures and emission measures as free parameters, temperature of the distribution, together with its uncertainty 

Results for the temperature of the coolest component {KT\) (standard deviation), that is 
were plotted versus time (Fig. |7^) and count-rate (Fig. [8^). 

No change in KT\ was noticed. Moreover, the normalized none , aai<;i \r 

fe ' KT 1 — 0.296 ± 0.016 keV. 

histogram of the time-weighted values of KT\ (Fig. [9p) 

showed an approximately Gaussian distribution peaked at the 

temperature of the coolest component of the 2-T model that This value can be considered the best estimation for the 

described the quiescent state. A Gaussian fit provided the mean temperature of the coolest component. Since KT 1 was constant 

during the observation, we fixed it and computed the fits 
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Table 3. Results derived from the time-resolved spectral analysis after fixing the hydrogen column density, the global coronal 
abundance, and the different temperatures (see text in Sect. 14.11 ). Columns show: (1) label referring to the time-region where 
the spectrum is included (see Figs. [TJ and [5]); (2) time-coverage; (3), (4) and (5) emission measures, from the coolest (EM{) to 
the hottest (EMj) thermal components; (6) reduced y 1 an d degrees of freedom; (7) and (8) X-ray flux and luminosity in the 
0.5- 10.0 keV band. 

(I) Region <2) Time slan - Time cnd (3) £Mi (4) £M 2 (5) £M 3 < 6 ^ M (d.oi) (7) /x <8) log L x (erg s~') 

(ks) (10 52 crrT 3 ) (10 52 crrT 3 ) (10 52 cm" 3 ) 

,+0.08 
-0.08 

+0.11 
-0.11 

+0.12 
-0.12 

+0.41 
-0.39 

+0.34 
-0.35 

+0.25 
-0.24 

+0.36 
-0.35 

+0.25 
-0.23 

+0.22 
-0.21 

+0.15 
-0.14 

+0.23 
-0.22 

+0.09 
-0.09 

+0.09 
-0.09 
+0.14 
-0.14 

+0.33 
-0.32 

+0.31 
-0.30 

+0.23 
-0.23 

+0.32 
-0.30 

+0.27 
-0.26 

+0.28 
-0.27 

+0.28 
-0.27 

+0.27 
-0.26 

+0.18 
-0.17 

+0.16 
-0.15 

+0.09 
-0.09 



L0 


0.0 


-3.6 


1 A-3+0-09 
1 - D -'-0.09 


i 4^+0.08 
^ -0.08 


0.10; 


FA 


3.6 


-5.6 


1.583H 


1 ^8+ 012 
1.38_ 12 


0.42; 


FB 


5.6 


-7.4 


1-67!^ 


1 54 +012 

i - JH -0.12 


0.42; 


FB 


7.4 


-7.7 


i Q9+0.38 
1 - 7z -0.37 


1 69 +0 37 

1.U7 3S 


1.17; 


FB 


7.7 


-8.2 


1 55 +0 - 32 

1 - JJ -0.31 


9 Qo+0.31 
^••"-0.30 


2.05; 


FB 


8.2 


-8.9 


2 00 +0 - 26 


9 i T+0.24 
^' -0.24 


0.92; 


FB 


8.9 


-9.3 


1 60 +0 - 33 

1-u -0.32 


1 76 +031 


1.65; 


FB 


9.3 


-9.9 


i o 1+0.26 
i - 81 -0.25 


i Q7+0.24 
L - y 1 -0.24 


0.50; 


FB 


9.9- 


- 10.7 


1 qs+0.23 
7 -0.22 


1 Q9+0.21 


0.70; 


FB 


10.7 


- 12.2 


i 71+O.I6 
1 - /1 -0.15 


1 64 +0 14 

1 - UH -0.14 


0.63; 


FB 


12.2 


- 12.9 


i qc+0.23 
7 -0.23 


1 61 +022 
'■"'-on 


0.54 


LI 


12.9 


- 17.7 


1 HQ +0.09 
i -'°-0.09 


1 67 +008 

i,u '-0.08 


0.69; 


L2 


17.7 


-23.7 


i o 1+0.08 

1 - 81 -o.oa 


r Q4+0.08 
L - y -0.08 


0.92; 


FC 


23.7 


-25.8 


1 89 +014 


1 98 +0 13 

1 - 7o -0.13 


0.55; 


FC 


25.8 


-26.2 


O 19+0.35 

-0.33 


9 99+0.31 
^■^ -0.32 


0.91; 


FC 


26.2 


-26.8 


9 i 9+0.30 
A - L -0.29 


2 07 +0 27 

-0.27 


1.95; 


FC 


26.8 


-27.9 


Lacs*] 


9 14+0.21 

^ -0.20 


1.44; 


FC 


27.9 


-28.4 


i qc+o.30 

1 - 7J -0.29 


9 i c+0.29 
z - 1J -0.28 


1.08; 


FC 


28.4 


-29.0 


9 qq+0.28 


1 Q 1+0.25 
1,?I -0.25 


1.00; 


FC 


29.0 


-29.6 


9 9 1 +0.28 
-0.27 


1 62 +025 

1-0 -0.25 


1.47; 


FC 


29.6 


-30.2 


1 89 +0 - 28 


1 98 +027 


0.99; 


FC 


30.2 


-30.8 


2 00 +0 - 27 


2 07 +0 - 26 

^■"'-0.25 


0.69; 


FC 


30.8 


-32.1 


9 1 c+0.19 
* ,J - J -0.19 


1 90 +0 18 


0.54; 


FC 


32.1 


-33.4 


2 00 +017 

^■""-0.17 


1 70+O.I6 

1. ' J_Q ] 5 


0.53; 


FC 


33.4 


-36.7 


1 Q9+0.10 


1 67 +009 

1 -"'-0.09 


0.12; 



(lO^'ergcrrT 2 s" 1 ) 


1.44(150) 


1.52 


29.38 


1.11 (209) 


1.77 


29.45 


0.97 (200) 


1.78 


29.45 


0.96 ( 85) 


2.39 


29.58 


0.87(161) 


3.17 


29.70 


0.97 ( 98) 


2.55 


29.61 


1.02 (111) 


2.61 


29.62 


0.89(114) 


2.14 


29.53 


1.21 (166) 


2.27 


29.56 


1.19(200) 


1.98 


29.50 


0.95 (147) 


1.99 


29.50 


1.22 (275) 


2.06 


29.51 


1.06 (264) 


2.37 


29.57 


1.36 (234) 


2.20 


29.54 


0.99(121) 


2.64 


29.62 


1.12(175) 


3.16 


29.70 


1.12 (184) 


2.81 


29.65 


0.91 (142) 


2.64 


29.62 


0.82(155) 


2.57 


29.61 


1.19(160) 


2.63 


29.62 


0.99 (156) 


2.47 


29.59 


0.96 (149) 


2.38 


29.58 


1.05 (148) 


2.24 


29.55 


1.15 (190) 


2.08 


29.52 


1.73 (139) 


1.74 


29.44 



again (second iteration) allowing variations in the other two 
temperatures and in the emission measure of all the three 
thermal components. No change was then detected in KT2 (see 
Figs-LTb and[8j)). This temperature also appeared to be uniform 
throughout the observation. Using a Gaussian distribution to 
fit the normalized histogram of the time-weighted values of 
this parameter (Fig.^), we obtained 

KT 2 = 0.88 ± 0.06 keV, 

which is clearly consistent with the temperature of the hottest 
component of the 2-T model that described the quiescent state. 
In addition to KT t , KT2 was also fixed in the model and the 
fits were repeated (third iteration) leaving free the temperature 
of the hottest component and the three emission measures. 
Figs.[7j; and[8j; show the values obtained for KTj, versus time 



and count-rate, respectively. It was not possible to constrain 
this temperature for the segments with an emission level equal 
or very similar to that of the quiescent state (first and last 
time-intervals). However, we cannot discard the existence of 
a little amount of hotter material during these time-intervals, 
with a contribution not strong enough to be detected. No clear 
pattern was detected in the time-evolution of KT3 during the 
observed flares (see Fig. [7};). In fact, the most surprising result 
is that no significant enhancement in this temperature was 
found during any phase of the observed flares when comparing 
them with non-flaring intervals (see also Fig. [8};). Furthermore, 
when accounting for the calculated uncertainties, all the values 
of KT3 are compatible with the average value obtained by 
fitting a Gaussian distribution to the normalized histogram of 
the time-weighted values of KT3 (Fig. [9}:), as occurred for the 
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Fig. 10. Time-evolution obtained for the emission measure of Fig. 11. Emission measure of the three thermal components vs. 
the three thermal components. Note that the time-resolution is count-rate, 
the same as that shown in Fig. [5] 

other two temperatures, that is measures as free parameters when doing the fits. Results of the 

best-fitting models are given in Table[3](note that the first and 

AT 3 = 1.87 + 0.30 keV. last time-intervals were also fitted with the same model). All 

the calculated fits can be considered good-quality models in 

Therefore, we fixed the three temperatures (KT\, AT 2 and terms of the^ ed test statistic. The time-evolution of the three 

KT3) and computed the last iteration leaving all the emission emission measures is plotted in Fig.[l0] Fig. [TTI complements 
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this figure giving the emission measure of each component ver- 
sus the count-rate. 

Our analysis shows that the occurrence of flares may be 
explained with the increasing of the emission measures related 
to the hottest temperatures, and in particular to the third one, 
while the dominant temperatures remain unchanged. 

4.2. Loop modeling 

Although stellar flares are not spatially resolved, we can 
infer the size of these structures by assuming that they are 
produced by the same basic physica l mechanisms as solar 
flares and using flare loop models (se e Kopp & Poletto, 1984; 



■ ■ ~. .... ~- ■ *~ ' ■ ' * . . L^^±m-^L^^LL2Z±±^D LJ^_^__Jl 

White et al., 1986; Poletto et al., 1988; van den Oord & Mewe, 

monlL„ii „ 1 „i I N rwi lo J Li I N r,r,iL Itt„...i„. ^ „i 



1989; 



1995 



Palla vicini et all Il990t ISerio et allll99lHHawlev et all 
iReale et all 119971: lor the review given bv Reale, 



2002. for an extensive discussion and comparison of these memo- 
Under the hypothesis of flares occurring inside closed coronal 
structures, and assuming that the heat pulse is released at the 
beginning of the flare, the decay time of the X-ray emission 
scales with the length of the loop which confines the flaring 
plasma. However, the presence of significant heating during 
the decay would slow down this phase and, therefore, the size 
of the flaring loop would be overestimated. It was shown that 
the slope (£) of the path of the flare decay in the density- 
temperature (log tip - log Tp) plane mainly depends on the 
heating decay tim e djakimiec et all 1 19921; ISylwester et al. , 
Il993t IReale et all Il997h . lRealeetal.1 dl997h derived an 
empirical formula using hydrodynamic simulations of single 
semi-circular flaring loops with constant cross-section, and 
including the effect of the heating in the decay: 



3.7 x 10-4F(<D 



F(0 > 1 



(2) 



where L is the loop half-length (in cm), r max the loop maxi- 
mum temperature (in K), tq the e-folding decay time derived 
from the light curve (in seconds), and F(0 a non-dimensional 
factor (larger than one) which accounts for the heating in the 
decay. The slope of the decay path in the density-temperature 
diagram is maximum (~ 2) if the heating is negligible dur- 
ing the decay, a nd minimum (~ 0.5) if the heating domi- 
nates this phase djakimiec et al. . 1992 ). Eq. (01 was success- 
fully tested on reso lved solar flares observed with Yohkoh/SXT 
(Real e et al. I ll997h and has been further extended and applied 



to several stellar flares observed with a wide variet 

J 



tectors dReale & Micelal Il998t iFavata et all l2000alblkl [200 It 
Maggio et all 12000b jGudel et all 1200 1[ IStelzer et all 2002[ 
Brig gs & Pvel \20oA (Reale et all |2004 IPiUitteri et all I2OO5I 
Fran ciosini et all 20061) . The correction factor F(£) needs to 
be calibrated for each detector, since it depends on the spec- 
tral bandpass and r esolution. Fo r EPIC-PN observations, the 
expression given by lReald d2007l) has to be used: 



of de- 



F(0 



+ q a 



(3) 



where c a - 0.51, ( a - 0.35, q a = 1.36. Eq. 0} can be used for 
slopes in the range 0.35 < f < 1.6. 



For stellar flares the density can rarely be directly mea- 
sured. Nevertheless, since reasonably the soft X-ray flare vol- 
ume is approximately constant during the flare development, 
the square root of the emission measure can be used as a proxy 
for the density (note that EMp - J n e riu dVp ~ n 2 Vp for a to- 
tally ionized hydrogen plasma). 

As described in § 14.11 the analysed flares (FB and FC) 
only involve the two hottest components of the 3-T model used 
to fit the spectra. We obtained the flare contribution of each 
one of these components to the total emission measure (.EM^f 
and EM^f) by simply subtracting from EM2 and EM3 their 
corresponding value measured for the quiescent state (L0). 
Therefore, the total emission measure of the flaring plasma 
(EM-p) is given by the sum of EM^y and EM^ p. A mean tem- 
perature of the flaring plasma (Tp) was also calculated for 
each spectrum as a weighted average of the temperature of 
Jhe two hottest components in the 3-T model, where the re- 
spective weights were EM2.p/EM F and EMt, p/EM f . The re- 
sulting log EMp 2 - log Tp diagram is plotted in Fig. Q~2] for 
the decay of FB and FC. A linear fit to the data provided a 
slope of f = 0.18 ± 0.27 for FB and f = 0.45 + 0.08 for FC, 
indicating the presence of sustained heating during the decay 
of both flares. Since the slope for FB is compatible with the 
lower asymptotic value for which Eq. <[3j can be applied, we 
obtained an upper limit for the loop half-length of this flare by 
using ( - 0.45. The maximum temperature in the flaring loop 
(T^max) can be derived from the maximum observed tempera- 
ture (r F max ). The observed temperature is a kind of average 
temperature of the flaring loop, which is therefore lower than 
the loop maximum temperature. The relationship between r max 
and Tp max depend s on the instrumental response. For EPIC-PN, 
Reale et all d2004 reported: 



= 0.1 84T, 



1.130 



(4) 



Thus, r raax resulted to be 32 X 10 6 K for FB, and 30 x 10 6 K 
for FC, with uncertainties of 16 %. 

Using the quantities given above, the values of Td shown 
in Table 12 and Eqs. d2j and (01, we obtained L < 7 x 10 9 cm 
for the half-length of the loop involved in the flare FB, and 
L ss (14 ± 8) x 10 9 cm for that involved in FC. 

5. Discussion and conclusions 

We have presented a detailed study of the X-ray variability that 
the binary star CC Eri showed during an XMM-Newton ob- 
servation. The great sensitivity of the EPIC-PN CCD camera 
allowed us to perform time-resolved spectroscopy of two flares 
weaker (flux increases of factors 1.5 — 1.9) than those typically 
analysed in the literature. The maximum luminosity of these 
flares in the 0.5 - 10.0 keV band was * 2.6 x 10 29 erg/s. 

Results show that, during the analysed observation, 
the corona of CC Eri is well-described by a 3-T model 
with constant temperatures (values of 3, 10 and 22 MK) 
and time-variable emission measures. The emission mea- 
sure of the coolest component (EM\) ranges between 
1.6 - 2.3 x 10 52 cm 4 , the one of the second component (EMt) 
between 1.4-2.3 x 10 52 cirT 3 , and that of the hottest com- 
ponent (EMt,) varies between 0.1 -2.1 x 10 52 cirT 3 . That is, 
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Fig. 12. Density-temperature (log «f - log T?) diagram of the variable component (see text) for the decay phase of flares FB (left 

1 /2 

panel) and FC (right panel). EM F ' has been used as a proxy for the density. Numbers indicate the time-evolution of the values. 
The point at the flare maximum (max) is also plotted. For each flare, the dashed line shows the best linear fit to the decay data. 
The derived slopes {() are also given. 



EM\ reaches 1.4 times the value measured for the quiescent 
state, EM2 varies up to a factor of 1 .6, and EMj increases up to 
21 times its quiescent value. Therefore, the hottest component 
dominates the X-ray variability observed on CC Eri (see also 
Fig. [Tol l. The changes in the emission measure of the coolest 
component could be due to the variable aspect of the binary 
during its orbital motion. However, the changes in the emis- 
sion measure of the two hottest components are clearly cor- 
related with the variations observed in the light curve during 
the detected flares (see Figs. [TOl and ITTb- The results obtained 
for the two hotte st components resemble those reported by 
Reale et al. d200lh for a sample of solar flares that covered a 
wide range of intensities and physical conditions. They found 
that the EM(T) distribution of all the flares in the sample fol- 
lows a common evolution path: it starts low but already at a 
relatively high temperature (~ 10 7 K independently of the flare 
intensity), it grows toward higher EM values during the rising 
phase and then decreases during the decay, maintaining always 
a more or less constant shape and peak temperature. In other 
words, they found that the height of the EM(T) distribution 
is clearly variable during the different phases of solar flares, 
while the width and peak temperature of the distribution suf- 
fer from much smaller changes. As far as we are concerned, 
this is the first time that no evident temperature variations are 
observed in a time-resolved study of stellar flares (see Figs. [7] 
and [8). Besides, the fact that the temperatures that character- 
ize the flaring plasma coincide with those of the quiescent state 
is consistent with a recurrent idea in the literature, that is: the 
quiescent emission of magnetically active stars may be pro- 
duced, to a large fraction, by continuous flaring activity (see 
§[TJand references therein). Thus, the light curves of multiple 
small flares overlap, and only larger flares stand out from the 
quiescent level. Besides, without temperature variations, differ- 
ences in the emission measure between time-intervals with di- 
verse count-rate, such as those observed for L0, LI and L2 (see 
Fig. [JJ and Table [3), can be identified with differences in the 
total coronal volume occupied by the flaring plasma. Density 



variations may also account for the changes in emission mea- 
sure. However, we expect that, if they were significant, they 
would be driven by local heating injections or leaks, and there- 
fore coupled to temperature variations but, as discussed above, 
temperature variations are not detected. 

The slope of the path in the density-temperature diagram 
indicates the presence of significant heating during the de- 
cay of the analysed flares, and leads to an upper limit of 7 
and 14 x 10 9 cm, respectively, for the half-lengt h of the flar- 



ing lo ops. From the orbital solution of CC Eri (lAmado et al 



2000) and the stellar radius of its components, we obtain a 
minimum distance between their surfaces of ~ 1 .4 x 10 11 cm. 
Even the largest loops involved in the detected flares only cover 
~ 10 % of this distance. Therefore, it is likely that the flaring 
loops are included all in the corona of one of the stars, and do 
not result from the magnetic fields bridging the two stars. In 
addition, bearing the spectral type of the stellar components of 
CC Eri in mind, the ratio between the half-length of the flar- 
ing loops and the stellar radius (/?») would be 0.1 and 0.3 if 
the flares were produced by the K7.5Ve star, or 0.2 and 0.4 
if they were produced by the M3.5Ve star. In both cases, the 
loop height above the stellar surface (2L/n for vertical loops 
with semi-circular geometry) is about 0.1 - 0.3 R*, implying a 
relatively compact flaring corona. These results are in agree- 
ment with t he loop half-leng t hs d erived for other Me dwarfs 



(L < 0.5 jL: lFavata& Micela, 2003). 



Assuming that at the flare peak the flaring loop is not 
far from a steady-state condition, and given that the derived 
loop half-lengths are significantly smaller than the pressure 
scale heighn , we ca n apply the so-called RTV scaling laws 
(Rosn er et all Il978h . These relationships link the pressure 



2 The pressure scale height is defined as h p = 2&r max /(jug), where 
yi is the average atomic weight and g is the surface gravity of the star. 



Therefore, h p as 5000T mllx 



Kg/ go)- Taking the spectral type of the stel- 
lar components of CC Eri into account, h p >7x 10 10 cm for both FB 
and FC. 
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Table 4. Parameters derived for flares FB and FC. 





Flare FB 


Flare FC 


_ T D (a) 


iOJU 


J7UU 


_^x,max (iv erg s ) 


Z.u 


z.o 


c Fv . , C10 33 erpt 


0.75 


1.5 


a 


0.18 ±0.27 


0.45 ± 0.08 


'J^ (10 6 K) 


32 ± 5 


30 ± 5 


£L (10 9 cm) 


<7 


14 + 8 


Ijp (10 3 dyn cm" 2 ) 


> 1.7 


0.74 ± 0.34 


*n e (10 11 cm- 3 ) 


> 2.0 


0.88 ± 0.43 


y F (10 30 cm 3 ) 


< 0.7 


3.2 


(erg s" 1 cm" 3 ) 


> 3.7 


-0.8 




~ 38 


~ 19 


iS(G) 


> 210 


> 140 



" e-folding decay time derived from the light curve. 

* Luminosity at the flare peak (0.5 - 10.0 keV band). 
Total radiated energy (0.5 - 10.0 keV band). 

d Slope of the decay path in the density-temperature diagram (see 
details in §|42}. 

e Maximum temperature in the loop at the flare peak (see § 14.2b . 
' Half-length of the flaring loops (see § 14.21 1. 

g Maximum pressure in the loop at the flare pe ak (estimated from 
the loop scaling laws given bv lRosner et al.Lfl978h . 

* Maximum electron density in the loop at the flare peak. We have 
assumed a totally ionized hydrogen plasma (i.e., p = 2n c AT raax ). 

' Volume of the flaring region (note that EM F « n\ V F for a totally 
ionized hydrogen plasma). 

' Heating rate per unit v olume at the flare pe ak (estimated from the 
loop scaling laws given bv lRosner etailll978h . 

* Number of loops needed to fill the flare volume. We have assumed 
a loop aspect r/L = 0.1 for a single loop. 

' Magnetic field (p < p B = B 2 /&7r). 



(p) and the heating rate per unit volume (En) with the loop 
half-length and the loop maximum temperature. Table [4] lists 
all these quantities for the observed flares. We have also esti- 
mated the electron density under the assumption of a totally 
ionized hydrogen plasma, obtaining n e ~ 10 11 cm" 3 . This is 
compatible with values expected for a plasma in coronal con- 
ditions. In order for the electron density to be consistent with 
the EMf measured at the peak of the analysed flares, the flar- 
ing volumes should be as large as 10 30 cm 3 . From the pressure 
of the flaring plasma (see Table [H, we infer that the minimum 
magnetic field required to confine the plasma at the flare peak 
is ~ 210 G for FB and ~ 140 G for FC. 

To satisfy the energy balance relation for the flaring re- 
gion as a whole, the maximum X-ray luminosity must be lower 
than the total input energy rate at the flare peak (H = EuVp). 
The rest of the input energy is used for thermal conduction, 
kinetic energy and radiation at lower frequencies. From the 
analysis of a large opt ical flare on the M3.5Ve star AD Leo, 
Houdebine et al.1 d 19931) concluded that the total kinetic energy 
during the event was of the same order as the radiated energy. 
For both FB and FC, we obtain that Lxmax is about 10 % of 



H, compatible with the X-ray radiation being only one of the 
energy loss terms during the detected flares. This value is in 
agreement with those reported for solar flares, where the soft 
X-ray radiation at the p eak only accoun t s for 10-20% of 
the total energy budget dWu et all Il986l) . iReale et alj d2004 
found a similar percentage (~ 15 %) for a flare observed on 
the M5.5Ve star Proxima Centauri, while Favata et"al] |2000c) 
estimated a higher value (~ 35 %) for an extreme X-ray flare 
detected on the M3.5Ve star EV Lac. 

If the detected flares were produced by a single loop, its as- 
pect (J3 = r/L, where r is the radius of the loop cross-section 
derived from the volume and loop length) should be 0.6 for 
FB and 0.4 for FC. Such a large cross-section is not observed 
on solar coronal loops, for which typical values of f3 are in 
the range 0.1-0.3. Therefore, we suggest a more realistic 
scenario consisting on flaring structures made up of several 
similar loops. Assuming f3 = 0.1, FB and FC occur in ar- 
cades composed of ~ 38 and 19 loops, respectively. Similar 
structures are also observed to flare on the Sun. For exam- 
ple, the Bastille day flare (2000 July 14) is an intense solar 
flare (GOES class X6) that occurred on a curved arcade with 



some 100 post-flare loo ps dAschwanden & Alexanderi 12001 



Reeves & Warrenl 12002). Stellar analogues have been proposed 
for the Me dwarf Proxi ma Centauri , where the best description 
for a flare analysed by IReale et al.1 d2004l) shows the presence 
of an arcade made up of ~ 5 loops; a nd also for the you nger 
(~ 100 Myr) G9 dwarf ZS 76, where IPillitteri et alj d2005l) es- 
timated 20 - 30 loops for the flaring arcades. Therefore, we 
conclude that events like solar arcade flares may be a common 
phenomenon on stars, in wide generality. 
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